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ABSTRACT 

A model is presented utilizing a Hamiltonian with equal spin singlet and triplet pairings based on quantum field 
theory and green function formalism, to show the correlation between the superconducting and ferromagnetic 
order parameters. The model exhibits a distinct possibility of the coexistence of superconductivity and 
ferromagnetism, which are two usually incompatible cooperative phenomena. The work is motivated by the 
recent experimental evidences of a long-range magnetic order below the superconducting phase temperature in 
a number of ternary or pseudo-ternary compounds of the lanthanides. The theoretical results are then applied to 
show the coexistence of superconductivity and ferromagnetism in the lanthanide compound HoMo 6 S$. 

KEYWORDS: Spin singlet and triplet state, Superconducting and ferromagnetic order parameters, Green 
function. 
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I. INTRODUCTION 

The study of magnetic superconductors has gained tremendous importance since the discovery of the 
ternary compounds of the lanthanides, in which a long-range magnetic order is observed below the 
superconducting phase temperature [1, 2]. In fact, ErRh 4 B 4 was the first ferromagnetic superconductor in which 
superconductivity was found to exist in a small temperature interval with modulated ferromagnetic phase as had 
been observed in a detailed study by Sinha et al. [3]. Later ferromagnetic superconductivity was experimentally 
confirmed in other compounds such as URhGe, UCoGe and ZrZn 2 [2, 4, 5]. 

The first observation of a zero resistance in the ferromagnetic state of HoMo 6 S8 was reported by Lynn 
et al.[6] and then by Genicon et al. [7]. Subsequently, neutron diffraction experiments carried out on singlet 
crystals of HoMo 6 S8 by Rossat et al. [8], further confirmed the above observations. They found that for T<0.60k, 
with slow cooling, a ferromagnetic phase appears and then the magnetic intensity increases and afterwards 
saturates at above T=0.4k, indicating that below this temperature, the thermal fluctuations are negligible. The 
material has a superconducting transition at T= 1.82k and a magnetic transition at T=0.67k in the vicinity of the 
reentry to the normal conducting state (Lynn et al.)[6]. 

The above exciting discovery stimulated a lot of interest in the study of coexistence of 
superconductivity and ferromagnetism. The two cooperative phenomena have been considered until very 
recently to be incompatible in view of the fact that conventional superconductivity favors spin singlet 

The problem of understanding the true nature and mechanism of superconductivity coexisting with 
ferromagnetism is so complex that a complete solution is still lacking [13, 14]. There has been suggestion that 
the coupling of longitudinal and transverse spin fluctuations plays an important role and should be included in 
the study of the problem [15]. Abrikosov [16] and Mineev et al. [17] have pointed out that S-wave 
superconductivity may result from the electron interaction mediated by ferromagnetically aligned localized 
moments. In the present paper, we start with a model Hamiltonian which incorporates not only terms of the BCS 
type but also an additional new term representing interaction between conduction electrons and localized 
electrons. The results of the calculations clearly show that over a certain temperature range superconductivity 
and ferromagnetism can co-exist. 



www.theijes.com The IJES Page 17 



Coexistence of Superconductivity and Ferromagnetism in... 



II. MODEL HAMILTONIAN OF THE SYSTEM 

The system under consideration consists of conduction electrons and localized electrons, between 
which exchange interaction exists. Within the framework of the BCS model, the Hamiltonian of the system can 
be written as: 

H=Hi+H 2 +H 3 (i) 

where, 

H i = S K/ff e K &w&w + Eu ej SIjB^ (2) 

represents, the single particle energies of the conduction and the localized electrons, measured relative 
to the chemical potential, d^^&n $■') and h- ii7 (fi^) are the creation (annihilation) operators for conduction and 
localized electrons respectively. 

*z = - J ^ fi kT <*! K j n _ K At ( 3 ) 

is the BCS type electron -electron pairing interaction term due to exchange of phonons. And, 

(h.c.= hermitian conjugate ) 

is the new term describing the interaction between conduction and localized electrons with the coupling constant 
u. the Hamiltonian in (1) will be used to determine the equations of motion in terms of the Green function. 

III. COUPLING OF SUPERCONDUCTING AND FERROMAGNETIC ORDER 

PARAMETERS 

The double-time temperature dependent retarded Green function is given by (Zubarev) [18]: 

G r (t-t") =«A(i);B(t")» 

= -i9it-t r ~)<[Ait),Bit")]> (5) 

Where A and B are Heisenberg operators and 8(t — t } is the Heaviside step function, defined as, 

{0 , if t < t 
which is related to the Dirac delta function as, — \.9{t — t )] = Sit — O (7) 

at 

Using the Heisenberg equation of motion (with h = 1), 
_dA(£) 



dt 



= [AttJM] 



along with the property (7), it is straightforward to obtain the following equation of motion foiG r it — O 

i— G r it - t) = Sit - O < [Ait) Jit}] > +« [A(t),H],B{t'}] » (8) 

at 

The Fourier transformation G r (,cJ) is given by 

G r (t -t")- I G r (tu}exp [-icait - 0]da (9) 

Taking the Fourier transform of (8), we get: 
uG r (v) =< [A(t}, § it"}] >^ +« [A(tl HI B(t-)] »^ (10) 

From (10), it follows that 

us « d kT fi + „ »= 1 +« [a^ Hi E+ a » (11] 

where the anti -commutation relation, 
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few **x*'}= *■»:'***'< ( 12 > 

has been used. To derive an expression for «. o kT , fi + H » , we have calculate the commutator [2^ ^ 
, using (2), (3), and (4). After some lengthy but straightforward calculations and using the identities 
[A, BC] = {A, B}C - B{A, C] md[AB,C] = A{£,C}-U,C}B " (13) 

we arrive at the following results: 

lz KV ff 3 ] = - £ % al^fi-pi^, (i4) 



[a^flj] = ^ "\,m«!*i *ff*i. 



J fRT 

Plugging (14) in to (1 1), we get 

td « a Kl , o 1 "^ »= 1 + e k « o kT ,o + kT » -2^ % «aI Ki a_ pi ap^ a 1 "^ » 
P 

k 

Using Wick's theorem for factorization ( Schwabl) [19], we can reduce (15) into the following form, 
k. ~e x > «**,**« »=1~ Cft - Y) « at Kl ,a^ rf » (16) 

where 

A= F^ « a_ p4 ,V »= ?£ « *V4 T » (17) 

y = u £ « b„, b^ » = u £ « £+ , S+, » (IS) 

ft ft 

are related to the superconducting and ferromagnetic order parameters respectively, and are assumed to be real. 
To determine <:< ajjy, o^f » in (16), we use (10), to obtain 

w « at u , r kT » =« [ ^.i*].^ » (i9) 

where (12) has been used, 

Proceeding in the same manner, we can reduce (19) into the form: 
fc. + ej « al^a*^ »= - (A - y) « ^.^ » (20) 

Eliminating <:< a,^, a^ » from (16) and (20), we obtain, 

- (A -y) 
« a+ .«* »= — — (21 1 

To take into account the temperature dependence of order parameters, we shall rewrite (17) and (18) as: 

K 

where 8 = 

For mathematical convenience, we replace the summation in (22) by integration. Thus 

T = J JV(0)de (24) 

where A r (0) is the density of states at the Fermi level. 
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Substituting (24) and (22) in (21), we obtain, 

For further simplification, we use the Matsubara frequency (Allen and Dynes 20 ), 

ftt -* itu n where, w n = ! — -^- (26) 

Using (26) and the fact that for effective pairing interactions, -ftwj < £ < hto b , we have from (25) , 

■^r ( A _ 7 ) 



A=2VN(0)sY I 



(2n + l) 2 ir 2 + /? 3 £ 3 



de, (27) 



where E 2 = f z + (A - y) 2 



Finally, using the Poisson summation 



And defining a = V T iV(0), we can rewrite (27) as 



=tanh (/? v 1 e 2 + (A - y)V2 J de :'2S : 

« J D V e + ( A - /) 

From (29), it clearly follows that the order parameters A and y , for superconductivity and magnetism are 
interdependent. 

We now consider the equations of motion for the Green function for the localized holmium f- 
electrons, which are responsible for magnetism. From (10), we can write, 

ta « b u , &+„ » = l+« [ b n .H].B i „ » (30) 

Substituting for H and following exactly the same procedure, we arrive at the relation, 



aZ-ef-Al 



where ^= £ ^ « 4 ^ » TO 

substituting (31) in (23), and using (26) and (28), we arrive at the following result 
y f to * 1 






(33) 



From (33), it is again evident that the order parameters A and y are interdependent, as was the case from (29). 

It is, therefore, possible that in some temperature interval, ferromagnetism and superconductivity 
can co-exist, although one phase has a tendency to suppress the critical temperature and the order parameter of 
the other phase. 
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IV. DEPENDENCE OF THE MAGNETIC ORDER PARAMETER ON THE 
TRANSITION TEMPERATURE FOR SUPERCONDUCTIVITY AND FERROMAGNETISM 

To study how y depends on the superconducting transition temperature T c , we consider the case, when 
T^O, orP^oo 



We can then replace 
tanh(£VE 2 + (A-y}V2] - 1 
in (29) and get, 
A _ r 1 ^* (A-y) 

« ~ h Je' + ta-ry 

or 



:d£ 



H-lh"^) C34) 

since \J~) 5®* 1, the above equation reduces to 

U(i-jW5!i) 

a v A J \fl — yt 

from which it immediately follows that, 

1 
&-Y = (2ftwj)exp(-— - — yi-r) (35) 

"U - 7aJ 

from the BCS theory, the order parameter A, at T=0 for a given superconductor with transition temperature T c is 
given by 

2A(C0 =3.5k s T c (36 

using this result in (35), we obtain 

/" 1 
y = l,7Sk B T c - 2frcj b gxp\ y (37] 



H (l 



l.75k B T c ' 



Once we know a, we can solve (37) numerically to draw the phase diagram for y and T c . 
To estimate a, we consider the case, 

T^T C 
which implies, A— >0 
From (29), we then have 

\ f f >^b J 



i i ■—■ i . , 

, tanhQ?Jc z + y a /2)d£ 



f^b y , 

-Jim tanh(gJc z + (d-r) z /21df 

*W D A ^i + (i - y )*} ^ v 

=; L -/ 2 (38) 

Putting £ : + y a = F : , 

we can write 



J 1= -tanhO?F/2)riF 

Using the result 

f tcmhx T T In j: 

j__ ^ =0nJc)(ianhx) | o -J— —i, 

D D 

and the fact that, at low temperature, taiih (- — — } -* 1 
we have from (39), 



(39) 
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^<m-o 



(40) 



where S is the Euler constant having the value S = 1.78 (Hsian) [21] 
we can write (40) as, 

^=111(1.14-— ) 

Using L' Hospital's rule, it is easy to show that 



-I 



(v-3)- 



2( e 3 +y 3 ) 

which can be neglected since y z is very small. 
Substituting (41) in (38), we then obtain 



dE 



a \ k E Tj 

This implies, 



1.14ft w 6 
T f = exp 



(-1) 



(41) 



(42) 



which can be used to estimate exp(— -1 for HoMo 6 S8, using the experimental value T c =1.82 k for this 

compound (Lynn et al.[6]). 

To study how y depends on the magnetic transition temperature T m , we consider (33). Neglecting A 3 , 
and proceeding as before, it is easy to show that, 



Y as -(aA)ln(1.14^-) 



which gives, 



From (36), we can estimate 4(0) « 0,2HmeV for HoMo 6 S8, using the known value of T c . 
we can use (43) to draw the phase diagram for y and T m . 



(43) 
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Fig. 1: Co-existence of superconductivity and ferromagnetism in HoMo 6 Ss 
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V. RESULTS AND DISCUSSION 

In Fig. 1, we have presented the theoretical curve of the magnetic order parameter y as a function of the 
superconducting temperature T c . For this purpose, we have used (eqn. 37) which has been numerically solved 
using the relevant parameters for HoMo 6 S8- In the same figure, we have also plotted the curve of y as a function 
of the magnetic transition temperature T m , using (eqn. 43). This curve is found to be almost linear up to the 
temperature T m = 0.67 k, which is the experimental value for HoMo 6 S8- 

From Fig. 1, we observe that T c decreases with increase in y, whereas T m increases with increase in y. 
The superconducting and ferromagnetic phases, therefore, oppose each other. However, the present model 
shows that there is a small region of temperature, where both the phases may appear together. This is indicated 
by (SC+FM) in the figure. 

We have thus demonstrated that a simple model based on a Hamiltonian which takes into account the 
spin interactions between conduction electrons and ferromagnetically ordered localized electrons, can explain 
the co-existence of singlet superconductivity and ferromagnetism in HoMo 6 S8- 
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